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£3 ' Abstract 

We consider a model with two parallel (positive tension) 3-branes separated 
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00 ■ by a distance L in 5-dimensional spacetime. If the interbrane space is anti- 

deSitter, or is not precisely anti-deSitter but contains no event horizons, the 
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effective 4-dimensional cosmological constant seen by observers on one of the 
branes (chosen to be the visible brane) becomes exponentially small as L 
grows large. 
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Recent observational data |J indicate that there is a positive cosmological constant in 
the universe, which, compared to the Planck or the electroweak scale, is many orders of 
magnitude smaller than expected within the context of ordinary gravity and quantum field 
theory. This is the well-known cosmological constant problem 0. Here we propose a solution 
to this problem in which the cosmological constant becomes exponentially small compared 
to all the other scales in the model. 

The scenario is a variant of the Randall-Sundrum model || . Let us consider two parallel 
3-branes (or two stacks of 3-branes) in 5-D spacetime separated by a distance L in the fifth 
dimension which is not necessarily compactified. Let the brane on the left be the visible 
brane, with positive brane tension cr . The brane on the right has positive brane tension 
d], / (To, generally. Assume that the bulks outside the branes are anti-deSitter(AdS) spaces, 
with cosmological constants —A; and — A r , respectively. There may be numerous solutions of 
the 5-D Einstein equations for the behavior of the bulk between the branes, but let us focus 
attention on the simplest, an AdS space with bulk cosmological constant — A m . (We discuss 
other possibilities briefly below, and in more detail elsewhere.) We shall see that A m may 
be expressed as a function of L, A;, A r , a , <Jl and the 5-D gravitational coupling constant 
k 2 = 8ttG. Thus, we regard A m as a derived parameter; more generally, in a dynamical 5-D 
spacetime, A m may evolve with time as branes move together or apart. For this simplest 
model for the interbrane spacetime, we can determine the effective 4-D cosmological constant 
A e ff in terms of A/, A r , <t , <Jl, L and k 2 . We show that for large L, 

A eff ^F(K,A h A r ,a ,a L )e- 2aoL (1) 

where F is independent of L and the positivity of the 4-D Newton constant Gjy requires 
a = K 2 a /3 — \J k 2 Ai/Q > 0. For a L ^> 1, the effective 4-D cosmological constant becomes 
exponentially small in this simple model without any fine tuning of the parameters. In a 
S 1 /Z2 orbifold version of the model, the only bulk is between the branes, and 

Aeff « 2a — — e . (2) 
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These results remain true even if ctq, which is simply the vacuum energy density of the 
4-D quantum field theory (which includes the standard model of strong and electroweak 
interactions) on the visible brane, changes due to phase transitions or other dynamics on 
the brane. We shall also comment briefly on solutions where the bulk between branes is 
not a pure AdS space. Demanding the absence of an event horizon between the branes, we 
again find that A e // to be vanishing small for large L. This suggests that an exponentially 
small effective cosmological constant may be a robust property of two-brane models with 
large interbrane separation. One is tempted to speculate that the pure AdS solution is a 
stable fixed point as the branes move apart and the 3-brane universe expands. 

When a L is not large, A e ff(L) is more complicated than eqs.(|l]) and (Q). Furthermore, 
at small separations, we expect additional nongravitational brane-brane interactions, but 
for large separations, it is reasonable to assume that the brane dynamics is dominated by 
pure gravity as described here. In a more realistic situation, the matter density on the 
visible brane (and dark matter on the other brane) should be included, and the separation 
distance L should be treated as a dynamical variable. For slow-moving branes, a Born- 
Oppenheimer-like approach is valid, and the spacetime evolves quasistatically from one 
nearly time-independent solution to another along the sequence of two-brane models. The 
result resembles the quintessence picture [|J and will be discussed elsewhere. Below, we 
derive eqs.(|l]) and @. 

Consider a pair of 3-branes at y = and y = L with different positive brane tensions 
<r and <j l . Let the bulk cosmological constant be —A; for y < 0, — A r for y > L. Choose 
the energy-momentum tensor for the bulk between the branes to have the diagonal form 
Tab = (— Ao, A, A, A, ip), where Xo(y), X(y) and if)(y) are functions of y. Choosing Ao = A 
allows us to use the metric ansatz 

ds 2 = dy 2 + A{y) [-dt 2 + e im 6 ij dx i dx i ] . (3) 

The G 5 component of the Einstein's equation Gab = k 2 Tab = SttGTab is satisfied trivially, 
while the Goo and the G55 components give, respectively, 
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AH ott2 o 2 

T = — + —i^(-y) + KvMvMl - y ) 

+A r e(y — L) — a 5(y) - a L 5(y - L)] 

Al\2 ATT2 2k 2 

-j) = — + —[Ai®(-y) + ^(yMyML-y) 

+A r Q(y-L)}. (4) 

The Gij component does not yield additional equations (alternatively, with the above metric, 
this component requires Ao = A). In the bulk between the branes, we have (differentiate 
if>(y) in eq. || and manipulate) 

2A' 

^' + ^-(^-A) = 0. (5) 

It is clear that if if) = A then X(y) — if)(y) = A m , independent of y. 

Let us first consider this particularly simple case, in which the space between the branes 
(0 < y < L) is AdS, with a cosmological constant —A m that will be determined from 



the solution of the 5-D Einstein equations. Define ki = Jk 2 Ai/Q, k m = JK 2 A m /6 and 



k r = ^K 2 A r /Q. Outside the branes, the solution for the bulks is 

A(y) = » 2 ^Ah(« + m)] {y < 0) 
A(y) = gWg Uy + y m )} {0<y<L) 

A(y) = ***?M±M {y > L) , (6) 

where yi, y m and y r are constants. These bulk solutions are similar to those in Ref ||. 
Continuity of the metric at the branes implies that 

smh 2 (kiyi) sinh 2 (/c m 7/ m ) 



h 2 P 

sinh 2 [fc m (L + y m )} _ sinh 2 [£y(L + y r )] 

U2 ~ h.2 

The jump conditions at the two branes are 

k\ cosh(kiyi) k m cosh(k m y m ) 



(7) 



sinh(%,) smh(k m y m ) 
k m cosh[A; m (L + y m )] _ k r cosh[k r (L + y r )} 
sinh[k m (L + y m )} sinh[k r (L + y r )} 
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go 

to- (8) 



where go = ft 2 ero/3 and qi = k 2 o"l/3. 

The expansion rate seen by observers on the brane at y — is H(0) = H/ JA(0), where 

H 2 _ _ 



A(0) smh 2 (kiyi) smh 2 (k m y m ) 

[kl ~ (k + qomi ~ (k ~ Qo) 2 ] 



(9) 



with A; 2 ^ — (ki ± go) 2 > or < 0, in agreement with Ref ||, which uses a slightly different 
approach. Similarly, the expansion rate seen by observers on the brane at y = L is = 



H/JA(L), where 



-ff 2 A; 2 fc 2 ^ 



A(L) sinh 2 [A; r (L + y r )\ sinh 2 [k m (L + y m )} 
[k 2 m -(kr + qL) 2 }[k 2 m -(k r -q L ) 2 } 

Aq\ ' UUj 

with fc 2 ^ — (k r ±q L ) 2 > or < 0. We can rescale t so that A(0) = 1, and the Hubble constants 



on the two branes are, respectively, H(0) = H and H(L) = H/yA(L). Note that although 
eqs.([|) and (|T0|) appear to determine the expansion rates on the two branes completely in 
terms of local quantities (i.e., the local brane tensions, and bulk cosmological constants just 
outside each brane), the values of these quantities on/near the two branes are connected via 
k m and y m . 

The 4-D Newton's constant Gn can be determined by introducing a small matter density 
p to the visible brane, that is, go — > Qo + n 2 p/3. Requiring the Hubble constant H to have 
the standard form 

rT n 87VGN / a \ /-, -,\ 

H 2 ~^-(Aeff + P+...) (11) 

yields 4nq GN = n 2 a k l Positivity of requires a = q — ki > 0; to be specific, let 

us consider 

< k m < g - k. (12) 

Since the expansion rates H and H(L) given in eqs.(|9]) and fllHD depend on k m , our goal 
is to express H and H(L) as functions of L and the parameters ki, k r , go and qi- This 
requires an expression relating L and k m \ from eqs.(^) and @ we find 



^ L = sinh " 1 (^)) +Sinh " 1 (^zy) 

= sinn 

V^-^ + *) 2 K-(^-<?o) 2 ] 
+ sinh- ( - 2kmqL ) (13) 

Here, eq.flllf) is regarded as a relation that determines k m in terms of L, q , qi, h and k r . 
In general, H(0) ^ H(L). Because of the condition (|12|), — > as A; m — > «o = Qo — h 
from below. This means that the expansion rate as seen by observers on the visible brane 
becomes exponentially small as L increases, 

H 2 ss 4a 2 e- 2{aoL - c) (14) 

where ao > 0, and sinh(C) = k m /H(L). As k m — > «o from below, -ff(-^) approaches a 
L-independent constant given by eq. (^0| ) and so does C. This implies that A e ff becomes 
exponentially small as L increases. 

In the symmetric case, where k\ = k r = k and q = qi = q (and let a = q — k = a ), we 
have H(L) = H(0) = H, and y m = —yi = y r = —L/2. The constant C in eq.flHD becomes 
C = aL/2, so 

H 2 « Aa 2 e~ aL + ^^p. (15) 
3q 

Thus, K e ff still decreases exponentially with L, but slower than in the nonsymmetric case. 

To this point, we have concentrated on spacetimes that are noncompact in y, but similar 
results can be derived for the compactified case. First, we may choose to identify ki = k r = 
k m and derive k m and H in terms of L. Next, we can compactify the y direction and further 
perform a Z 2 orbifold, with one brane sitting at each of the two fixed points (y = 0, L). This 
S l /Z 2 orbifold model is particularly simple, since there is only one bulk space between the 
branes sitting at the two end points. This is an expanding (non-supersymmetric) version 
of the Horava-Witten model ||10|| , with branes with tension o"o = 3go/^ 2 at y = and 
o"l = 3qi/K 2 at y = L, separated by AdS space with bulk cosmological constant —6k 2 /n 2 . 
The solution is 
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H 2 

A(y) = —smh 2 [k(y-sy )} (16) 

where s = +1 for y > and s = — 1 for y < 0. Because of the symmetry of the model, we 
only need to consider the jump conditions at y = and y = L, which are 2k j qo = tanhfcyo 
and 2k /qi = tanh[/c(L — yo)], respectively. Combining the jump conditions implies 

g _ tanhfcL - q L /2k 

2k ~ 1 - (q L /2k)tanhkL' ^ ' 

if qL/2k = ±1, then go/2& = =Fl irrespective of fcL, but for 7^ 1, go/2& -> 1 as 

kL —> 00. According to our viewpoint, eq.flPTD determines k given q , q L and L. 



The expansion rates on the branes are H(0) = H/ y A(0) and H(L) = Hj JA(L), where 
H 2 (0) = k 2 [(q /2k) 2 - 1] and H 2 (L) = k 2 [(q L /2k) 2 - 1]. For large values of q L we find that 

h 2 (o) « go(gx + go) e _ go x (lg) 

Thus, for q^L ^> 1 and qL > Qo, the cosmological constant on the y = brane is exponentially 
small. Moreover, although eq.([0|) may appear singular as qL qo, in fact if 2 (0) ~ q 2 e~ goL ^ 2 
in that case. If qL = —qo, as considered in Ref ||, then y = ±00, so H(0) = H(L) = 
and therefore |go| — 2fc for any finite non-zero L. Our model is qualitatively different, since 
it involves two positive tension branes, and interprets the bulk cosmological constant as a 
parameter derived from the brane tensions and separation. 

Notice that A(y) given by eq.(^[) may vanish at certain values of y. We view points where 
A(y) = as event horizons since it takes a test particle leaving either brane an infinite time 
to reach them according to observers on the brane from which it is launched. (There are no 
spacetime singularities at these points, since all derivatives of A{y) are finite.) We do not 
worry about event horizons outside the branes (y < and y > L), which recede to large 
distances from the branes as L increases, but an event horizon between the branes might 
be worrisome. To avoid an event horizon between branes, we need y m < L + y m < 0, which 
can be true if k 2 — q\ — «g > 0. This condition can be satisfied when the two branes and 
the outside bulks are not identical, but in the symmetric case, avoiding a horizon between 



the branes requires a < 0, which violates the positivity of Gat. So, in the symmetric case, 
we will have an event horizon right between the branes even when they are close together 
(L — ► 0). (Non-gravitational interbrane interaction could alter this conclusion for very small 
L, but not in general.) In the orbifold case, although we can allow arbitrary values of go 
and qi, horizons are only avoided in the model if L < y Q , which implies q^ < 0. Fortunately, 
horizons can be avoided altogether if the space between the branes is not pure AdS. 

There is nothing sacred about AdS space between the branes (see |ll|]), and we might 
expect non-AdS behavior to be the rule rather than the exception, particularly at different 
stages in the evolution of a slowly- changing 5-D spacetime. Generally speaking, we can find 
infinitely many families of solutions with different X(y) and ip(y) in the regions between 
branes. As a first stab at a more complicated two-brane model than we have presented here, 
we have studied the case where A = A m is a constant, but 4>(y) = A m — K/A 2 (y), which 
satisfies eq.(|5[) for y— independent K. In this case, we do not get a relation between k m and 
L (because ip(y) ^ A m ). The question of whether H 2 is small is related to horizon formation 
between the branes. If horizons are absent, then H 2 must be exponentially small for large L. 
If a horizon may only appear between branes at large L, then if 2 must also be exponentially 
small. In a dynamical solution with separating branes, we might imagine that H 2 — > as 
t — > oo, and no horizon forms. 



To be specific, let us consider the orbifold model in some detail. Here, we find (for K > 
and A(0) = 1) 



with s chosen as for AdS. In this case the jump conditions at y = and y = L are 




(19) 




qh_ 
Ik 



( 



H 2 \ sinh[2A;(L - y )] _ 
2P/ A(L) cosh2ky ' 



(20) 



using the first of these relations in the second implies 



S 



H 2 _ [(go + Ql)] cosh2A;L - (1 + qoq L /2k) sinh2£;L 
~k ~ sinh2A;L - (q L /2k)(l - cosh2fcL) 

+ (21) 



q L ~2k 

Horizons are guaranteed to be absent in this model if 

H 2 1 



< 1 ~ 2e~ 2hyo (22) 

2k 2 cosh2/b/ -l ' 



Using eq. (pll) in the jump condition at y = implies e 2kyo ~ e fe y (7l(<?o — 2k)/q (qL — 2k), 
and so horizons are absent if 

2k q L -2k K J 

in the general case, qi ^ qo- If ?o — then e~ 2kyo ~ e _fei and horizons are avoided if 
g /2A; - 1 < 2e- kL . In either case, the expansion rate is exponentially small if there are no 
horizons in the solution. 

The dynamical treatment of such a separation modulus may follow that discussed in 
Ref [F2|] and is under investigation. For large separations, it is reasonable to assume that 



the brane-brane interaction is dominated by pure gravity as described above. At small 
separations, we expect brane-brane interactions beyond pure gravity. In Ref [Q] where an 
AdS-Schwarzschild solution is considered, the deviation from the pure AdS space red-shifts 
away rapidly as the universe expands. This leads one to conjecture that the pure AdS 
solution between the branes is a stable fixed point. If so, as branes move apart and the 
universe expands, the bulk solution remains horizon-free, and approaches pure AdS space 
asymptotically. Clearly this issue needs careful investigation. 

One may even envision a multi-brane scenario, where the separation distances between 
branes play the roles of various scalar fields, as suggested by string/M theory. In partic- 
ular, the separation of two nearby branes may play the role of an inflaton [|13|], while the 
separation distance between two far-apart branes may play the role of a scalar field simi- 
lar to quintessence p|. Immediately after inflation, the separation of two nearby branes is 



stabilized by short-range brane-brane interactions [fl2| , [l3[1 , when the effective cosmological 



constant is small compared to the radiation density. Then a third brane starts to move away 
while the 3-brane universe is expanding, yielding an exponentially small A e jj. 

We thank Eanna Flanagan and Horace Stoica for useful discussions. This research is 
partially supported by NSF (S.-H.H.T.) and NASA (I.W.). 
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